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Abstract. The 7 Cephei planetary system is one of the most interesting systems due to several 
reasons: 1 .) it is the first planet candidate detected by precise radial velocity (RV) measurements that 
was discussed in the literature (Campbell et al. 1988); 2.) it is a tight binary system with a 20 AU; 
and 3.) the planet host star is an evolved K-type star. In Hatzes et al. (2003) we confirmed the 
presence of the planetary companion with a minimum mass of 1 .7 Mj up at 2 AU. In this paper we 
present additional eight years of precise RV data from the Harlan J. Smith 2.7 m Telescope and 
its Tull Coude spectrograph at McDonald Observatory. The 900 d signal, that is interpreted as the 
presence of the giant planetary companion, is strongly confirmed by adding the new data. We present 
an updated orbital solution for the planet, which shows that the planet is slightly more massive and 
the orbit more circular than previous results have suggested. An intensive high-cadence week of RV 
observations in 2007 revealed that 7 Cep A is a multi -periodic pulsator. We discuss this issue within 
the context of searching for additional planets in this system. 
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INTRODUCTION 

7 Cep is a nearby (d = 13.8 pc), bright (V = 3.2) Kl IV-III star, that has the distinction 
to be the very first candidate for an extrasolar planetary companion detected by precise 
RV measurements that was discussed in the literature. The pioneering RV survey at 
the CFHT carried out by Bruce Campbell and Gordon Walker discovered that 7 Cep 
is not only a single-lined binary, but also that it shows periodic low-amplitude RV 
variations that might indicate the presence of a giant planetary companion with a period 
around 900 days (Campbell et al. 1988). However, the same team later cast doubt on 
the existence of this planet (Walker et al. 1992). They argued that this signal could be 
caused by intrinsic variability of the star itself. More recently, Hatzes et al. (2003) used 
over 14 years of precise RV data obtained by the McDonald Observatory planet search 
to demonstrate that the 900 d signal stayed coherent in phase and amplitude of the entire 
time span. We also showed that there is no significant variability in the chromospheric 
emission level and line bisectors based on the McDonald spectra of 7 Cep. We concluded 
that the 900 d signal is indeed due to a planetary companion with a minimum mass of 
1.7 Mj up orbiting at a = 2.1 AU. 

The detection of this planet in a relatively tight binary system with a ~ 20 AU lead 
to a host of follow-up investigations ranging from dynamical studies (e.g. Dvorak et 
al. 2003, Haghighipour 2006), to planet formation modeling (e.g. Kley & Nelson 2008, 
Jang-Condell et al. 2008), to additional work combining existing RV data with astro- 
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FIGURE 1. Almost 30 years of precise RV data for y Cep A. The four different data sets are from the 
CFHT (triangles) and from the 3 different phases of the McDonald Observatory planet search (circles, 
diamonds, filled circles). The solid line displays the RV orbit due to the stellar secondary. 



metric measurements (Torres 2007), and to direct imaging campaigns (Neuhauser et 
al. 2007). In particular the combination of the last two papers significantly improved 
the parameters of the 7 Cep stellar binary system. The orbital period of the binary is 
67 ± 1.4 years with an eccentricity of 0.41 and a semi-major axis of 20.2 ±0.7 AU. 
The primary has a mass of Ma = 1 .4 ± 0. 12 M and the secondary is an M dwarf with 
M B = 0.41 ±0.02 M . 

In this paper we present our results based on additional eight years of precise RV 
measurements from the planet search program at the Harlan J. Smith 2.7 m telescope 
(HJST) and its Tull Coude spectrograph at McDonald Observatory. 

UPDATED ORBITAL AND PLANETARY PARAMETERS 

After the announcement of the planet in 2003, we have continued to observe 7 Cep as 
part of the long-term precise Doppler survey at the HJST. We are currently in phase III 
of the program that utilizes the entire available spectral bandwidth of the iodine cell that 
we use to measure precise differential RVs. Phase I and II used either different reference 
lines (e.g. telluric lines) or instrumental configurations. A detailed description of the 
different phases of our program can be found in Hatzes et al. (2003). Figure 1 displays 
the entire RV data starting with the CFHT results and ending with the McD phase III. 
These data cover the periastron passage of the binary orbit (solid line). With these new 
data in hand, we can confirm or falsify the planetary hypothesis by testing if the residual 
signal stayed constant in period, phase and amplitude. 
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FIGURE 2. Lomb-Scargle periodogram of the RV residuals after subtracting the binary orbit (upper 
panel) and the window function of our observations (lower panel). The peak at 900 days has an estimated 
false-alarm-probability of 10~ 30 . 



Fig. 2 shows the Lomb-Scargle periodogram of the RV residuals after subtracting the 
large binary motion. The 900 d signal is clearly present in our data and its statistical 
significance is higher than in 2003. The false-alarm-probability of the peak is « 10~ 30 . 
Our phase III spectra contain the Ca II H&K lines that can be used as proxies for 
magnetic activity of the star. We determine a Ca line S-index that corresponds to the 
level of chromospheric emission in the Ca H & K line cores. If any correlation between 
the S-indices and the RV data is found it would cast serious doubt on the existence of a 
companion. Fig. 3 contains the periodogram of our S-index measurements. The vertical 
dashed lines indicates the 900 d period for the planet. It shows that there is no isolated 
peak at the orbital period, further strengthening the case for the planet. However, there 
is a peak of moderate power at a significantly shorter period of 800 d. Although 
not highly significant, this peak could indicate the rotational period of the star. Our re- 
analysis of the CFHT FWHM measurements in Hatzes et al. also revealed a shorter 
period that can be the stellar rotational period. 

We conclude that the addition of the new data strongly confirms the presence of the 
900 d signal residual signal and that there appears to be no connection between this 
signal and intrinsic variability of the host star. 

We derived a simultaneous Keplerian orbital solution for both the binary and the 
planetary orbit using Gaussfit (Jefferys et al. 1988). The RV offsets between the four 
individual data sets were included as fit parameters. The only parameter that we held 
fixed is the period of the binary. We used 24650 days, the current best value following 
Torres (2007) and Neuhauser et al. (2007). Letting this parameter vary freely lead to 
shorter binary periods, which are inconsistent with the global solution based on historic 




FIGURE 3. Lomb-Scargle periodogram of the Ca II H& K line indices from our McDonald spectra 
(upper panel) and their window function (lower panel). The vertical dashed line displays the period of the 
planetary companion. There is some power in a neighboring peak at w 800 days that might indicate the 
rotational period of y Cep A. 



low precision RVs and imaging results. Despite that the 30 years of precise RVs cover 
the periastron passage and almost half of the entire orbit, it is not sufficient to derive 
the correct binary period using these data alone. The remaining orbital parameters 
for the binary are in good agreement with the published values: K= 1932 ± 3 ms -1 , 
e = 0.41 ±0.001, (O = 162±0.3 and T Per = 2448502± 8.5 days. The updated parameters 
for the planet are given in Table 1. The biggest difference to the previously published 
values are a higher minimum mass of 1 .85 ± 0. 16 Mj up (previously 1 .7 ±0.4 Mj up ) and a 
lower eccentricity of 0.049 ±0.034 (previously 0.12 ±0.05). Fig. 4 shows the 4 different 
RV data sets (after subtracting the binary motion) along with the planetary orbit phased 
to the 900 d period. 



TABLE 1. New orbital and planetary pa- 



rameters for 7 Cep Ab 


parameter 


value 


error 


units 


Period 


903.3 


1.5 


[days] 


K 


31.1 


0.97 


[ms- 1 ] 


e 


0.049 


0.034 




Tper 


2453227 


87 


[days] 


(O 


94.6 


34.6 


[deg] 


m sin i 


1.85 


0.16 


[M Jup ] 
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2.05 


0.06 


[AU] 



in 

> 
T3 



> 
OS 

T3 



100 
50 



- 



-50 



100 



1 — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — p 



"CFHT 




J i i i I i i i I i i i I i i i I i i i L 



1 — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — r 



"McD I 




J i i i I i i i I i i i I i i i I i i i L 




FIGURE 4. The four RV data sets phased to the orbital period (903 d) of the planetary companion. 
The residuals scatter around this new orbital solution are: CFHT: 15.4 ms~', McD I: 19.5 ms -1 , McD II: 
17.7 ms- 1 , and McD 111:7.7 ms- 1 



Searching For Additional Planets 

The residual RV scatter of « 8ms _1 is more than twice as large as the internal 
measurement uncertainties of 3ms 1 for the McDonald phase III data. This excess 
scatter could be a sign of additional, previously undetected planets in the system. Based 
on dynamical simulations, Haghighipour (2006) showed that low-mass planets with 
a < 0.8 AU should remain stable in this system. However, a period search in the residuals 
to our binary plus single planet fit did not (yet) reveal any significant signals (see Fig. 5). 

In September 2007 we used seven consecutive nights at the HJST to perform an 
intensive asteroseismological study of y Cep A. We collected over 1200 individual RV 
measurements and discovered that the star is a low-amplitude multi-periodic pulsator. 
The results of this observing run are described in more detail in Endl et al. (2009). We 
identified 12 different frequencies with amplitudes ranging from 1.5 to 4.1ms _1 . A 
section of one night from this campaign is displayed in Fig. 6. The oscillation spectrum 
of y Cep A has a large spacing of 15/1 Hz, which is proportional to the mass and 
radius of the star. Using a mass of 1.4 M we derive a radius of 4.9 R . This value 
is right in between the two interferometrically measured radius values of 4.79 ±0.06 R© 
(Nordgren et al. 1999) and 5.01 ±0.05 R (Baines et al. 2009). 

The oscillations of 7 Cep A are the cause of the observed excess scatter. Although we 
have currently no evidence for additional planets in the system, such a search is limited 
by the short-term variability due to the pulsations. To account for this additional RV 
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FIGURE 5. Lomb-Scargle periodogram of the RV residuals after subtracting the binary and the single 
planet model. All peaks have a false-alarm-probability (FAP) of > 1%. Therefore, we have no evidence 
for additional RV signals in our data. 




signal in our data, we have increased the error bars on our phase III data artificially by 
a factor of two before performing the orbital fits presented in this paper. Otherwise, the 
small formal errors of the phase III data, which do not take the oscillations into account, 
would have severely underestimated the true uncertainties of these data. 

The pulsation signal of a few mr 1 is present in the entire RV data set and it limits 
our ability to search for lower-mass planets that would produce RV signals with similar 
amplitudes. The exposures times for a bright star like J Cep are short enough so that 
we sampled these pulsations at random phases, rather than average them out. Therefore, 
as the next step, we plan to define and subtract the oscillation signal from the entire 30 
years of RV data. This should allow us to reach the highest sensitivity for additional 
planets in the 7 Cep system. 
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